Zero field spin polarization in a 2D paramagnetic resonant tunneling diode 
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We study I-V characteristics of an all-II-VI semiconductor resonant tunneling diode with dilute 
magnetic impurities in the quantum well layer. Bound magnetic polaron states form in the vicinity of 
potential fluctuations at the well interface while tunneling electrons traverse these interface quantum 
dots. The resulting microscopic magnetic order lifts the degeneracy of the resonant tunneling states. 
Although there is no macroscopic magnetization, the resulting resonant tunneling current is highly 
spin polarized at zero magnetic field due to the zero field splitting. Detailed modeling demonstrates 
that the local spin polarization efficiency exceeds 90% without an external magnetic field. 



The implementation of device components based on 
resonant tunneling diodes (RTDs) is one route towards 
the elaboration of a full semiconductor spintronics based 
technology scheme. While a ferromagnet /tunnel barrier 
spin injector [UJ3] produces a fixed spin polarization for 
each given magnetization state, dilute magnetic semi- 
conductors (DMS) can be used in II- VI semiconductor 
RTDs to implement spin selective tunneling at different 
bias voltages [3] . A caveat to this approach has been the 
paramagnetic nature of bulk (Zn,Mn)Se, which makes 
the application of an external magnetic field necessary 
for spin filter operation. This can be overcome by using 
the 0D states of self assembled quantum dots embedded 
in a DMS host material, since the microscopic magnetic 
environment of a dot allows for the formation of bound 
magnetic polaron (BMP) like states which lift the spin 
degeneracy for the tunneling electrons [S] . Such self 
assembled quantum dot structures have a rich resonance 
spectrum which typically occur over a broad range of 
bias voltages, limiting the controllability of device char- 
acteristics. Here we show that similar zero field splitting 
can be achieved in the much more reliable quantum well 
geometry. 

We investigate an all-II-VI RTD grown on a GaAs 
substrate. The active RTD region contains a 9 nm 
Zno.96Mno.04Se quantum well layer sandwiched between 
two 5 nm Zno.7Beo.3Se tunnel barriers. Proper contact 
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FIG. 1. Schematic conduction band profile of the resonant 
tunneling diode at zero bias, with the spin degeneracy in the 
quantum well lifted in the quantum well. 



layers are applied on each side of this structure to allow 
for measurements of transport through the layer stack 
(fig. [I]). The quantum well layer is made from a DMS 
that exhibits giant Zeeman splitting in an external mag- 
netic field, which is described by a modified Brillouin 
function [3J [7] with a pair breaking contribution at high 
magnetic fields J>L Lifting the degeneracy of the quan- 
tum well spin states with an external magnetic field al- 
lows the RTD to be used as a voltage controlled spin 
filter [3]. The I-V characteristic shows current peaks at 
two different bias voltages as long as the splitting is large 
enough to resolve the separate spin up and down reso- 
nances. 

The black lines in fig. [2] show I-V characteristics for 
measurements at 1.3 K from to 14 T. Similar results 
for fields up to 6 T have previously been successfully 
described [3] using a model based on taking the conduc- 
tance of a single spin channel to be one half of the B=0 
T curve, applying Brillouin splitting to the quantum well 
levels and recombining the contribution of the two spin 
channels into a total TV curve by using Kirchhoff 's laws. 
Such a model implicitly assumes spin degeneracy at B=0 
T, and obviously breaks down if that condition is not 
fulfilled. The data presented here, which include higher 
magnetic fields than available previously, suggest that 
a modified picture of the zero-field tunneling process is 
necessary. 

As shown by the blue lines in fig. [2] the data is sugges- 
tive of the peak splitting not vanishing at B=0 T. More 
importantly, the peak in the zero field I-V characteris- 
tic is also less symmetric than each of the split peaks 
at high magnetic fields, and the resonance in the zero 
field curve is much broader than that of the individual 
resonances in the 14 T curve. Both the asymmetry and 
the increased width of the peak in the B=0 T curve may 
be a consequence of this peak actually being comprised 
of two resonances occurring at somewhat different bias 
voltages. These considerations indicate the need for a 
different modeling scheme. 

Bottom-up approaches to modeling such data have 
been reported O [10], but these typically treat an ide- 
alized system ignoring considerations such as contact re- 
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FIG. 2. Fits (red dashed curves) to the I-V characteristics 
(black) at 1.3 K and at magnetic fields from to 14 T, applied 
perpendicular to the layer stack. Each curve is offset by 14 
fiA on the current axis for clarity. The blue lines are a guide 
to the eye to emphasize the apparent peak splitting at B=0 
T. 

sistances. In fig. [3] we show that such considerations 
are important. Fig [3^, gives the equivalent circuit of our 
real device in the two channel model, and includes mag- 
netic field dependent contact resistances i?f'^, an inter- 
face scattering term R SC at and a non resonant contribu- 
tion to the tunneling current, R og . The active region 
of the RTD is represented by the two diodes, one for 
each spin channel, and each with a voltage and mag- 
netic field dependent resistance Rb g in parallel. While 
the diode carries the resonant part of the current includ- 
ing the LO-phonon replica contributions, the background 
current through R og accounts for electrons tunneling off 
resonance through the double barrier region. 

To obtain an expression for the highly non-linear re- 
sistance resulting from the resonant tunneling transport, 
one normally assumes a Lorentzian shaped transmission 
at the resonance condition. The peaks in our experi- 
ment show a Gaussian line shape with a bandwidth much 
broader than the expected injector Fermi energy (mul- 
tiplied by the lever arm). Thus an additional broaden- 
ing mechanism, probably stemming from imperfect inter- 
faces at the active RTD region, dominates the resonance 
width 11 . Local potential fluctuations, caused by well 
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FIG. 3. a) Resistor model for the two spin RTDs in parallel, 
b) Measurement and fit for the 14 T I-V characteristic. The 
blue and red curves represent the current carried by the spin 
up and down species, adding up to the purple curve which 
is the fit to the measurement (black curve), c) Plot of the 
potentials and V* at each spin diode as a function of the 
applied bias voltage. 



width fluctuations |12j or inhomogeneous alloy or doping 
concentrations |13) . impose an additional in-plane con- 
finement for tunneling electrons thus creating OD type 
tunneling states, so called interface quantum dots [14j . 
Since our device is 100 2 /im 2 and these fluctuations are 
typically on a nm scale, we sample over an ensemble 
of these states in our vertical transport measurements. 
One can view this configuration as a large number of OD 
resonant tunneling diodes in parallel, each with its own 
resonance conditions. This results in a broadened Gaus- 
sian line shape [T2| for the overall resonant conductance 
feature. The LO-phonon replica are described by addi- 
tionally broadened Gaussian conductance peaks with re- 
duced amplitudes and an energetic separation from their 
respective spin-split resonance peaks of 31.7 meV, the LO 
phonon energy of bulk ZnSe [15] . 

For the non-resonant background current we use a 
transfer matrix model consisting of two tilted barriers 
where the quantum well was omitted in order to re- 
move resonant contributions. The potential drop over the 
quantum well region, which effectively lowers the second 
barrier, plays an important role, and is explicitly taken 
into account. The resulting transmission is proportional 
to the non-resonant tunneling of emitter electrons and 
fits well to the measurement at high bias voltage, where 
the contribution of resonant tunneling is small. 

Due to the contact resistances R\'^, the two voltage 
nodes Vf and Vj, in fig. [3^i are not necessarily at equipo- 
tential for a given applied bias voltage V app . Fig. [SJ; 
shows the potential at the points and Vj, as a function 
of V app - When a resonance condition is reached for ei- 
ther of the spin diodes, the resistance of that spin diode 
drops and the potentials across each of the diodes is al- 
tered accordingly. While we have experimental access to 
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V app , the transport theory for resonant tunneling only 
describes the active region of the device. Thus consid- 
ering the contact resistances is vital for fitting any RTD 
model to actual experiments. As an example, the result- 
ing fits for a magnetic field of 14 T are presented in fig. 

where contributions from both the spin up and down 
channels are shown as well as how they add up to produce 
a fit (purple curve) to the observed measurement (black 
curve). While the conductance of a resonant channel is 
perfectly symmetric on an energy scale, fig. (3]d+c show 
how in a real device, the contact resistances influence the 
shape of the resulting l-V app characteristics. 

As is clear from the circuit diagram in fig. [3^,, the total 
current traversing the device is given by 



I{V ap p) = ^ (^(yt) + 4 o( yt) + ^(yt)) + 

V±(aHV±) + at (V±) + ai g (V±)) (1) 



with a t - l (V 1 ~' 1 ) ocp T4 • 



exp 



(2) 



where er^, crj^ and er^Q are the conductances for the 
spin channels, the background contributions and the LO- 
phonon replica peaks respectively. I is the lever arm link- 
ing the energy scale in the quantum well to the diode bias 
voltages V^'^, Eq'^ is the energy between the spin levels 
and the conduction band edge, p^^ are fitting parameters 
for the amplitudes of the spin conductances (and thus 
yield the spin polarization) and are the variances 
of the Gaussians describing the energy level distribution 

for the spin channels. Equation 2 is also used for u\^ 

^ i — . ^ 
but with different variances ^£q, amplitudes p£q and 

4xo = 4 l + 31-7 meV. 

Our detailed model therefore consists of solving the 
equivalent circuit of fig. [3^, for an RTD with a spin split 
resonance and the associate LO-phonon replica. Since 
the zero field I-V characteristic is a superposition of two 
strongly overlapping peaks, the best starting point for 
the fits is the high magnetic field data, where one eas- 
ily can find the proper variances and of the 
resonant peaks and LO-phonon replicas. Starting at 14 
T, the TV characteristic for each magnetic field is fitted 
by adjusting p^, p^ and Eq. We also allow for a 
magneto-resistance effect in the contacts RJ'^ and in the 
scattering channel R sca t- By including a magnetic field 
dependence of T^q , we account for the small, experimen- 
tally observed field dependent broadening of the replica 
peaks. 

The resulting fits are shown as red dashed lines on top 
of the TV characteristics in fig. [2] while in fig. [3]d the 
contributions of spin-up and spin-down electrons to the 
14 T I-V characteristic are illustrated. One would a pri- 
ori expect a Brillouin function to describe the magnetic 
field dependence of the splitting [4 . The measurements 



shown in fig. [2] exhibit a very different behavior. At 
low magnetic fields we observe that instead of a spin de- 
generacy, the l-Vapp characteristic is properly fit only by 
allowing for finite splitting even at zero magnetic field 
|llj . We have previously observed such a remanent zero 
field splitting in the zero dimensional resonant tunneling 
states of self assembled CdSe quantum dots [H] . Here the 
quantum well is nominally a two dimensional object. As 
previously discussed, however, various inhomogeneities 
cause the current transport to be effectively mediated by 
a large ensemble of parallel paths each flowing in a lo- 
cal environment. The relatively low number of magnetic 
atoms influenced by each of these regions means that 
each will statistically have, on average, a net magnetiza- 
tion at zero field [5]. This effect is further enhanced by 
the presence of the spin of the tunneling electron [TS] . 

The energy separation between spin-up and spin-down 
peaks is 15 meV at B=0 T as determined by the fit. This 
energy is not necessarily the same as the splitting of the 
two spin states. As the measurement is always referred to 
the conduction band of the emitter, this energy difference 
is influenced by the different bias conditions needed to 
align each spin state to the emitter. 

The first peak at lower bias voltage is suppressed while 
the second peak is enhanced in the B=0 T l-V app char- 
acteristic. For a small energetic splitting in the reso- 
nant state, one would expect similar conductances for 
the two transport channels. A change in the confinement 
caused by the splitting will influence the amount of leak- 
age of the quantum well wave function into the emitter, 
while for each V app the resulting change in symmetry of 
the double barrier will affect the transmission [17]. A 
higher bias voltage will also drive more current at the 
same conductance. From transfer matrix calculations for 
the transmission probabilities of the double barrier we 
conclude that different biasing conditions alone cannot 
explain the magnitude of the effect on the amplitudes of 
the B=0 T spin currents. 

While the peak positions stay constant at intermediate 
fields because the Brillouin function saturates, above 8 T 
there is a clear reduction in the splitting of the peaks on 
the bias voltage axis. A reason for this reduction is likely 
the Zeeman splitting of the emitter electrons, since both 
ZnSe and (Zn,Mn)Se have a positive g factor and the 
resulting splitting AV res on the voltage axis is given by 
Al/,. es = I (dQW — 9e) ^bB, where I is the lever arm of 
the device and g<Qw and 9e are the effective g-f actors of 
the ZnSe emitter and the quantum well electron states re- 
spectively. From the fits we obtain a slope of 0.47 meV/T 
(-0.18 meV/T) for the spin-up (down) peak. The corre- 
sponding g- factors are g^ = 16.3 and g^ = 7.3, far greater 
than the bulk ZnSe value of 1.1 15j. Possible explana- 
tions for this increase in magnetic response include that 
tunneling electrons at the interface to the barrier can- 
not be treated in the free electron picture of a parabolic 
s-type conduction band, that there is a dilute Mn con- 
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FIG. 4. Current spin polarization as a function of applied 
bias voltage and magnetic field. 



centration in the emitter due to diffusion during growth, 
or that the energetic distance to the resonant quantum 
well state is altered by spin selective band bending of the 
emitter. The peak amplitudes are also strongly magnetic 
field dependent. The asymmetry of the effective g fac- 
tors for the emitter polarization suggests an effect that is 
linked to the resonance bias conditions. The two peaks 
occur at different bias voltages and therefore have differ- 
ent conduction band bending conditions. This bending 
changes the number of available electronic states for res- 
onant tunneling, strongly influencing tunneling currents. 
This factor can easily surpass the effect of Zeeman split- 
ting. The resulting effective g-factors are therefore not 
purely a result of the electron spin interacting with the 
magnetic field but also of the feedback mechanisms in- 
duced by changes in the potential landscape QT]. Differ- 
ent transmittances of the spin channels may also result in 
spin sensitive charge build-up in front of the barrier that 
can influence not only the amount of available states in 
front of the barrier, but also the bias voltage needed to 
attain the resonance conditions [18 . 

The amplitudes we obtain from the above fitting 
process give quantitative results for the spin polarized 
currents. To use this device as a detector for the emit- 
ter spin polarization one would need to link the emitter 
polarization to the amplitude of the traversing spin cur- 
rents. The increased and asymmetric magnetic feedback 
that is evidenced by the movement of the spin peaks sug- 
gests that other effects in addition to pure Zeeman split- 
ting of the emitter are involved [TT]. Therefore usage as 
a detector for the emitter spin polarization is difficult. 
Within our model it is possible to evaluate currents for 
spin-up and down electrons separately, thus allowing for 
a quantitative analysis of the polarization of the current 
traversing the device. Fig [4^,) shows the current spin 
polarization as a function of magnetic field and bias volt- 
age. Blue (red) indicates a spin up (down) polarization 



of the current. The non resonant background current is 
not spin selective and therefore the current polarization 
P c plotted in fig [4] is given by 



Pr = 



(3) 



The splitting of the spin levels in the external magnetic 
field and the changes in the amplitudes of the resonant 
peaks lead to a polarization above 90% for the spin-up 
peak at B=14 T, while the spin-down peak polarization 
decreases to below 60%. While a high degree of polariza- 
tion of both spin types can be achieved at all measured 
magnetic fields, counterintuitively, despite the paramag- 
netic nature of bulk (Zn,Mn)Se, the maximum polariza- 
tion efficiency is achieved without applying an external 
magnetic field, where 80% for spin-up and 90% for spin- 
down is observed, as evidenced by the FV curve of the 
two channels for B=0 T presented in fig. |4Jd). Similar- 
results for a second device and at various temperatures 
are presented in pTj . 

In summary, we have shown high spin polarizations can 
be achieved due to formation of BMP like states in the 
active RTD region. The resulting microscopic magnetiza- 
tion for the tunneling electrons lifts the spin degeneracy 
and provides two separate transport channels. Feedback 
mechanisms stemming from the influence of different bi- 
asing conditions both increase the energy splitting of the 
peaks and influence their amplitudes, resulting in high 
degrees of current spin polarization. Our model allows 
for good fits to the device characteristics and thus quan- 
titative analysis of the polarization. Not only does this 
model confirm the findings of reference |4j that the device 
can work as a voltage controlled spin filter at moderate 
magnetic fields, but it also establishes that the local spin 
polarization efficiency not only remains, but is even en- 
hanced in the absence of a magnetic field. 

We thank A. Slobodskyy and T. Slobodskyy for sample 
preparation and measurements, and the DFG SPP 1285 
Semiconductor Spintronics for financial support. 



[1] 

[2] 

[3] 

[4] 
[5] 
[6] 

m 

[8] 
[9] 
[10] 



S. A. Crooker et al, Science 309, 2191 (2005). 

R. Jansen and B. C. Min, Phys. Rev. Lett. 99, 246604 

(2007). 

R. S. Patel, S. P. Dash, M. P. de Jong, and R. Jansen, J. 
Appl. Phys. 106, 016107 (2009). 

A. Slobodskyy el al, Phys. Rev. Lett. 90, 246601 (2003). 
T. Dietl and J. Spalek, Phys. Rev. Lett. 48, 355 (1982). 
C. Gould el al, Phys. Rev. Lett. 97, 017202 (2006). 
A. Twardowski, M. Vonortenberg, M. Demianiuk, and R. 
Pauthenet, Solid State Comm. 51, 849 (1984). 
Y. Shapira el al, Phys. Rev. B 30, 4021 (1984). 
P. Havu el al, Phys. Rev. B 71, 235301 (2005). 
A. Saffarzadeh, M. Bahar, and M. Banihasan, Physica E 
27, 462 (2005). 



5 



[11] See attached supplementary data and analysis. 
[12] A. Zrenner et al, Phys. Rev. Lett. 72, 3382 (1994). 
[13] O. Makarovsky et al, Phys. Rev. B 81, 035323 (2010). 
[14] D. Gammon et al, Phys. Rev. Lett. 76, 3005 (1996). 
[15] Landolt-Bornstein, New Series, Group III, Vol. 22a, 

Springer-Series (1999). 
[16] T. Dietl and J. Spaek, Phys. Rev. B 28, 1548 (1983). 
[17] S. Luryi, Superlattice Microstr. 5, 375 (1989). 
[18] Y. W. Choi and C. R. Wie, J. Appl. Phys. 71, 1853 

(1992). 



6 



supplementary data 




FIG. 5. I-V characteristics (black lines) and fits (red dots) 
of a sample with 8% Mn and 6% thinner layers at magnetic 
fields from to 14 T (in steps of 0.5 T) and T=1.3 K. Each 
curve is offset by 14 fiA on the current axis for clarity. 

Fig. [5] shows a comparison of our model to experimen- 
tal data for a second type of sample. This device has 8% 
Mn instead of 4%, and all layers are 6% thinner that the 
device in the paper. The graph presents I-V character- 
istics at 1.3 K for magnetic fields from to 14 T (black 
lines). The fits (red dots) again agree well with the exper- 
iment, even though, due the high Mn content, the spin 
down resonance is merged with the spin up replica peak. 
All fitting parameters show the same magnetic field de- 
pendence as in the manuscript. As in the text, the blue 
lines are a guide to the eye. 

Fig. [6] gives the details of the fit to the zero magnetic 
field I-V characteristic, again showing zero field splitting 
with a spin up peak of reduced amplitude and a pro- 
nounced spin down peak. In this case, a spin polarization 
of 80% (92%) for spin up (down) electrons is achieved. 
We now consider the effect of temperature dependence, 
by analyzing data taken at 6 T for various temperatures 
ranging from 45 mK to 15 K. For each temperature Temp 
one can solve the equation 

Brillouin[6 T,Temp] = Brillouin[B e// ,1.3 K] 

to determine at which magnetic field R e ff a curve from 
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FIG. 6. Also the control sample is highly suggestive of zero 
field splitting. The I-V characteristic (black line) is fitted 
(purple line) by adding the spin up (blue line) and down (red 
line) currents. 



the 1.3 K dataset has the same level splitting in the quan- 
tum well as the 6 T curve at the given temperature. Fig. 
[7] presents the level positions of the resonant spin states 
for both the 1.3 K data set of fig. [5]and this temperature 
dependent measurement. The open symbols are for the 
1.3 K dataset, and the x-axis is then directly the magnetic 
field at which the measurement was performed. The solid 
symbols are for the temperature dependent data, plotted 
against R e ff &s described above. This comparison con- 
firms that the movement of the peak position is a result 
of changes in the band diagram, and not a result of any 
deformation of peak shape, as these would not be stable 
under the different environmental conditions. 



> 0.065 
•SL \ 

O 0.060 



g_ 0.055 
ID 

oJ 0.050 p 



, o0 ^ 0o0 ogOoo»oO»oo^ 00 o °oo OQ 



8 10 12 14 



FIG. 7. Energetic distance to the conduction band of spin 
up (circles) and down (triangles) levels in the quantum well. 
Filled symbols show results from fits at various temperatures 
at B=6 T, while empty symbols show fits at various magnetic 
fields at T=1.3 K. 



A similar analysis also is conducted for the amplitudes 
of the spin channel conductances in fig. [8j The polariza- 
tion values acquired from fits of measurements at various 
temperatures differ from those acquired from the mag- 
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FIG. 8. Change in amplitudes of the spin up (circles) and 
down (triangles) levels with applied external magnetic field 
(normalized to the B=0 T amplitudes). Filled symbols show 
results from fits at various temperatures at B=6 T, while 
empty symbols show fits at various magnetic fields at T=1.3 
K. 
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FIG. 9. Self consistent conduction band profile of the resonant 
tunneling diode at the approximate resonance biases. The 
inset shows the increase of tunneling states due to the applied 
bias voltage. 



netic field dependent measurements at T=1.3 K. After 
correcting for the level splitting using B e y/, the only vis- 
ible difference between the two sets of measurements are 
peak amplitudes. This is because a constant quantum 
well splitting is maintained in the two configurations, 
which then only differ in emitter polarization due to the 
applied external magnetic field. Fig. [8] shows that in- 
deed, the only field where both polarization values are 
identical, is at B e ff= 6 T, where both temperatures are 



the same. For B e // < 6 T the temperature measure- 
ments show higher polarization since, while the splitting 
in the quantum well is maintained constant, the external 
magnetic fields effect on polarizing the emitter produces 
a higher polarization of the spin current. The opposite 
is true for B e // > 6 T. The results of fig. [8]thus suggest 
that both the splitting in the quantum well and the po- 
larization of the emitter influence the spin polarization 
of the resonant current. 

As ment ioned in the text, the reason why a change in 
140 
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FIG. 10. Comparison of I-V characteristics at T=40 mK, 1.3 
K and 15 K. The resonance does not sharpen at low temper- 
atures. 



quantum well splitting already results in a change of cur- 
rent polarization is the change in the bias voltage needed 
to reach resonant conditions for each of the spin channels. 
Fig. [9] shows self-consistent calculations of the conduc- 
tion band profile at the resonance condition for the 4% 
Mn sample. As the inset shows, the Fermi energies differ 
by approximately 20% for the two resonance conditions. 
Since the maximum current flows when the quantum well 
level is aligned with the conduction band edge and is 
proportional to the cross-sectional plane A — irkp of the 
emitter Fermi sphere at constant E z this would result in 
a w44% change of the peak amplitudes. 

Lastly, in fig. [I0]we plot the B=0 T FV characteristics 
of 40 mK (red dots), 1.3 K (gray line) and 15 K (blue 
dots) , showing that in the absence of magnetic field, tem- 
perature does not have any influence on the I-V charac- 
teristic. This indicates that a much stronger broadening 
mechanism is at work in the device; namely the potential 
fluctuations at the quantum well interface. 



